Thermal symmetry of the Markovian master equation 
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The quantum Markovian master equation of the reduced dynamics of a harmonic oscillator coupled 
to a thermal reservoir is shown to possess thermal symmetry. This symmetry is revealed by a 
Bogoliubov transformation that can be represented by a hyperbolic rotation acting on the Liouville 
space of the reduced dynamics. The Liouville space is obtained as an extension of the Hilbert 
space through the introduction of tilde variables used in the thermofield dynamics formalism. The 
angle of rotation depends on the temperature of the reservoir, as well as the value of Planck's 
constant. This symmetry relates the thermal states of the system at any two temperatures. This 
includes absolute zero, at which purely quantum effects are revealed. The Caldeira-Leggett equation 
and the classical Fokker-Planck equation also possess thermal symmetry. We compare the thermal 
symmetry obtained from the Bogoliubov transformation in related fields and discuss the effects of 
the symmetry on the shape of a Gaussian wave packet. 
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I. INTRODUCTION 

In the dynamical formulation of a quantum system in 
the Hilbert space TL, a continuous symmetry is repre- 
sented by a unitary transformation U that commutes 
with the Hamiltonian H [1]. When we extend the dy- 
namics to the Liouville space C ~ Ti. x Ti, the symme- 
try of H is translated into a symmetry of the Liouville 
operator, defined by L = —i{H x 1 — 1 x H) for a sys- 
tem with unitary time evolution. However, the reverse 
is not true. The Liouville operator may contain a sym- 
metry that has no counterpart in the Hamiltonian [2]. 
For example, consider a simple harmonic oscillator with 
Hamiltonian Hq = wpa^a. This system possesses a non- 
degenerate discrete energy spectrum, £'„. Nevertheless, 
the energy spectrum of the corresponding Liouville op- 
erator Lq = —i{HQ X 1 — 1 X Hq) depends on energy 
differences Em — En , and therefore degeneracies may oc- 
cur [3]. The occurrence of degeneracies in the spectrum 
of the Liouville operator is a signature for the presence 
of a symmetry in C that has no counterpart in Ti. 

It often occurs that the system we are interested in 
is embedded in a larger system, for instance a thermal 
reservoir. The subsystem is an open interacting quan- 
tum system with a non-unitary (semigroup) time evolu- 
tion. We obtain the reduced dynamics of the system by 
tracing out the degrees of freedom of the reservoir. When 
memory effects are further neglected, we obtain the quan- 
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turn Markovian master equation (MME) [4]. This for- 
malism provides a framework for describing dissipative 
processes in various fields, such as quantum optics [5, 6], 
quantum Brownian motion [7, 8], quantum decoherence 
[9] and, more recently, quantum information theory [10]. 
The generator of the time evolution of the MME is a 
non-Hermitian collision operator with complex eigenval- 
ues [11, 12]. For different thermal states of the system, 
degeneracies may occur in the eigenvalues of the collision 
(or Liouville) operator. These degeneracies indicate the 
presence of a thermal symmetry that is the main subject 
of this paper. In contrast to the discussion in Ref. [2], 
in which the authors considered an exact symmetry S of 
the Liouville operator in the sense of [S, L] = 0, we con- 
sider a symmetry U of the collision operator K(b) up to 
a similarity transformation, K{b') — UK{b)W, where b 
labels the different thermal states of the system. 

A natural setting for formulating our discussion is pro- 
vided by the thermofield dynamics (TFD) formalism [13]. 
In this formalism, the extension to the Liouville space is 
achieved by introducing the operator A ^ a x 1 and the 
tilde conjugate operator ~ 1 x a for the annihila- 
tion operator a of the harmonic oscillator. A Bogoliubov 
transformation, which can be represented by a hyperbolic 
rotation, can then be defined to act on the pair of opera- 
tors (j4, A''') and its Hermitian conjugate pair [14]. Ther- 
mal effects of the reservoir enter our analysis through a 
temperature dependent angle in the hyperbolic rotation. 
As a result, the rotation connects different thermal states 
of the system; in particular, each temperature is related 
to the zero temperature that contains purely quantum 
effects (without thermal effects). 

A thermal symmetry may also exist in a classical sys- 
tem, though it occurs in a different form. This is true 
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for the classical Fokker-Planck equation [15] where the 
thermal symnietry results from a scale transformation 
on the phase space variables. We find that the Caldeira- 
Leggett equation [16] possesses a thermal symmetry as 
well. However, when we consider the thermal symmetry 
of MMEs derived under more general conditions [17-19], 
the thermal symmetry is no longer present. 

Since the Bogoliubov transformation preserves the 
canonical commutation relations between the position 
and momentum operators, we can introduce a symplec- 
tic structure into the theory. The thermal symmetry can 
then be formulated as an element of the Sp(2n,R) group 
acting on the phase space operators, which is often em- 
ployed in the studies of "Gaussons," or Gaussian pure 
states [20]. We also discuss how the generator of the 
thermal symmetry is related to the Lie algebra of the 
SU(1, 1) group utilized in quantmn optics [21] and com- 
pare the thermal symmetry to other canonical transfor- 
mation discussed in the literature [22] . 

Since the thermal symmetry connects different ther- 
mal states of the system, it also dictates the change in 
shape of the density function in the coordinate space as 
a function of the temperature. It is known that a den- 
sity function with Gaussian profile experiences simulta- 
neous stretch and contraction along the diagonal and off- 
diagonal directions separately in the coordinate space as 
a result of a change in temperature [23] . We will show the 
connection between these phenomena through the ther- 
mal symmetry. We conclude with a discussion of our 
results. 



II. THERMAL SYMMETRY IN MARKOVIAN 
MASTER EQUATION (MME) 

Consider a harmonic oscillator interacting with a ther- 
mal reservoir. When memory effects can be neglected, 
e.g. in the weak coupling limit, the time evolution of 
the reduced density operator p of the harmonic oscilla- 
tor is governed by a Markovian master equation (MME), 
dp/dt = Kp. The time evolution operator is [5, 24] 

Kp={Ko + Kd)p, (1) 

where 

Kop = -iuji^[a^ a, p] (2) 
KdP = ci{2apa) — a^ap — pa) a) 

-\- C2{2a) pa — aa) p — paa)) , (3) 

in which we have made the rotating- wave approximation 
for the interaction [25] and a and satisfy the usual 
commutation relation. All the dissipative effects of the 
dynamics are contained in Kj,. This MME generates a 
completely positive dynamical semigroup on the reduced 
dynamics [4]. 

In Eq. (2), is the natural frequency of the harmonic 



oscillator and ci, C2 in Eq. (3) are coefficients given by 




In the weak coupling approximation, the damping con- 
stant 7 is given by [24] 

7=/ duj2-K\^\v{uj)\^5{uj-uJo) = 2n\^\v{uJo)\^ , (5) 

where A denotes the coupling constant of the interac- 
tion and v{u}) is a model dependent form factor. The 
delta function in Eq. (5) arises from the contribution of 
the resonance pole at uJq. The thermal properties of the 
reservoir are encoded in the parameter b. In the weak 
coupling limit, b can be determined from the following 
expression 

/>oo 

76= / (hj2'K\^\v{ij)\'^b{uj)5{uj - u)q) 
Jo 

= 2TrX'\v{ujo)fb{iOo), (6) 
where we obtain 

b = b{<^o) , (7) 

after using Eq. (5). Assuming that the reservoir is in 
thermal equilibrium, the parameter b{u>) is given by 

bi^) = 1 + ^mJ^^ = I c°th (i/lu;/?) , (8) 

where /3 = l/fcsT, ks being Boltzmann's constant and 
T being the temperature of the reservoir. The physical 
range of b that corresponds toO<T<oois5<6<oo. 

We will now formulate our discussion in terms of the 
Liouvillc space by following the superoperator formula- 
tion in TED [13]. We introduce the (super-)operators 

A = axl, A"! = a'< xl. (9) 

and their tilde conjugate operators [32] 

A=lxa\ A^ = lxa. (10) 

These operators act on the density operator from the left 
according to, e.g., AAp = A{pa'^) = apa\ etc. Hermitian 
and tilde conjugations are idempotent operations [33], 

(At)t = A, (5)= A, (11) 

and they act on the product of two arbitrary operators 

X, Y as 

(7/Xy)t = ry*Ftxt , {riXY) = r}*XY , (12) 

where 77 is a complex number and * denotes complex 
conjugate. The commutation relations of the operators 
are 

[A,A^]^l, [I,At] = l, (13) 
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while all other commutators vanish. We also have the 
useful relation 



[A-A^,{A-A^)^]=0. 



(14) 



In terms of these operators, the uncoupled and dissipa- 
tive components of the collision (super-)operator become 

Ko = -icJn {A'' A - A'' A) 



.Wo 



{A + A^)\a - At) + {A + At) _ 



(15) 



and [34] 
a 2 



(^A-A^y{A + A^) - (A - At) (A + At) 
-57(A-At)^(A-At), 



(16) 



respectively. 

Now consider the unitary operator 



u{e) 

with the Hermitian generator 

G = i(AA- AtAt) 



(17) 



(18) 



Under the action of U by conjugation, the A operators 
are transformed into 

A' = UAU'^ = cosh OA - sinh 61 At , (19a) 
A't = C/Atj/t = - sinh OA + cosh 6»At , (19b) 

with similar expressions for their Hermitian conjugate 
partners. Eq. (19) is a Bogoliubov transformation which 
preserves the commutation relations in Eq. (13). This 
transformation can be represented as a hyperbolic rota- 
tion on the operator pair (A, At) by [14, 21] 



where the rotation matrix is 

m- 



A 
At 



cosh 6 sinh Q 
sinh 9 cosh 9 



(20) 



(21) 



The transformation on (At , A) is obtained by taking the 
Hermitian conjugation of Eqs. (19) or (20). It should 
be stressed that R{6) (or U) mixes the operators A and 
its tilde conjugate At, and hence it is a rotation in the 
Liouville space that has no counterpart in the Hilbert 
space. 

The effect of U on the operators Kq and can now 
be deduced by using the following linear combinations of 
Eqs. (19) 



A'TA't = e±^(ATAt). 



(22) 



The uncoupled component of the dynamics Kq has exact 
symmetry under [/, 

i^^(A,At) = UKo{A,A})U^ = Ko{A',A'^) 

= K^{Aj^), (23) 

that is, \U, Kq] = 0. This symmetry in Kq leads to the 
degeneracies in the spectrum of the operator Lq (= ^^^o) 
mentioned in the introduction. For the; dissipative com- 
ponent of the dynamics /ST^, the first line on the rhs of 
Eq. (16) has exact symmetry under U, but the second 
line, which is proportional to h, transforms up to a simi- 
larity transformation (or remains form invariant). Over- 
all, we have 



K'M. A^-M= KM', A'^;b) = KdiA, At ; 6') 



(24) 



where the transformed parameter b' is related to the angle 
of the rotation by 



e = ^Hb'/b). 



(25) 



For a reference thermal configuration b, the physical val- 
ues of 9 lie within the range — (ln2fe)/2 < 9 < oo, for 
^ < 6' < oo. As a consequence of Eqs. (23) and (24), 
the hyperbolic rotation connects the full collision opera- 
tors between reservoirs of different temperatures up to a 
similarity transformation 



K{b') = U{9)K{b)UH0)- 



(26) 



We note that in Ref. [2], the term "Liouville symmetry" 
refers to the existence of a symmetry S that commutes 
with the Liouville operator L (or the collision operator 
K), i.e. [L,S] = 0. For our system, while the thermal 
symmetry U{6) does not commute with K, it preserves 
the form of K up to a similarity transformation between 
different thermal states of the system (26). 

A quantum statistical system is described by a den- 
sity operator p defined on the Hilbert space Ti. As we 
formulate the dynamics in the Liouville space C, the sys- 
tem is now described by a density state vector jp)), which 
corresponds to p [3, 21]. We have used the double ket no- 
tation to emphasize the fact that the state vector lies in 
C. Under the thermal symmetry, the density state then 
transforms as 



\p'{b))) = U\p{b))) = \p{b'))) 



(27) 



and the time evolution is now governed by d\p{b')))/dt = 
K{b'Mb'))). 

The existence of the thermal symmetry is related to the 
degeneracies in the eigenvalues of the collision operator 
C{b) = iK{b) in the following way. Consider the non- 
Hermitian eigenvalue problem of C{b) [11, 12, 27], 

C{b)\ptnm = ^t.n\ptnm , (28) 

which has complex eigenvalues denoted by 

z^^ = ±rvjJo — i{m — n/2)^, m>n, (29) 
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where m,n = 0,1,2,3,..., and the decay rate of the 
(m,n) mode is (m — 71/2)7. Under the action of U, 
Eq. (28) is written 



nin 1 1 run 



ih') 



(30) 



Since the complex eigenvalue is independent of the 
thermal parameter b, degeneracy in can occur for 
different thermal states of the system. This degeneracy 
in the eigenvalues of the collision operator i^T is a conse- 
quence of the thermal symmetry. 

The transformation U can be interpreted as generat- 
ing a continuous change in the value of T (or h). For a 
transformation of the system from a lower temperature 
to a higher one, we have b' > b (equivalently, we can 
view this as a decrease in the value of H) with a posi- 
tive angle, 9 > 0. The system thus moves away from 
the quantum limit, while a transformation in the oppo- 
site direction with ^ < drives the system closer to the 
quantum limit. In the following discussion, we shall fix 
the value of h and interpret the transformation as gener- 
ating a change in the temperature. 

Under the U transformation, a system that is initially 
characterized by b is mapped into a different thermal 
state characterized by b'. Hence the term thermal sym- 
metry. Moreover, all finite temperature states are related 
to the T = (6 = i) state through the transformation. 
Since the T = state is a purely quantum state, this im- 
plies that for systems with a thermal symmetry, such as 
Eq. (1), thermal effects can be inferred from the effect of 
vacuum fluctuation by means of the thermal symmetry. 



where p = —id/dx is the dimensionless momentum. 
With similar expressions for the tilde operators, we find 
that the operators and their tilde counterparts separately 
satisfy the canonical commutation relations, 



[x,p\ =1. 



[X,P\ =1; 



(32) 



while all other commutators vanish. 
With the definition 



p{x, x) = {{x; x\p)) = {x\p\x) , (33) 
the A operators are represented by [12, 24] 



'{x; x\A\p)) = ^(x+-^] p{x, x) 



V2 



dx 



{{x;x\A^\p)) = (2^- ^) ^(^'^) 



(34a) 
(34b) 



in the coordinate basis, with similar expressions for the 
tilde operators. 

Prom Eq. (34), the pair of coordinates {x, x) transforms 
under the thermal symmetry just like the pair of {A, A^) 
operators, 



X J ^ ' \ X 



We also have the transformation equations 



(35) 



III. THERMAL SYMMETRY IN OTHER MMES 

In the previous section, the thermal symmetry is pre- 
sented in the superoperator form in Eq. (20), which 
clearly illustrates that the symmetry cannot be reduced 
to a symmetry in the Hilbert space. Other aspects of the 
symmetry can be discussed more conveniently in terms 
of the coordinate basis representation. For instance, the 
existence of a thermal symmetry in the classical Fokker- 
Planck equation becomes obvious when we compare it 
with the quantum MME in the Wigner's representation. 
The coordinate basis also proves to be more convenient 
when we extend our discussion to include other MMEs, 
where they acquire more simple forms. 

We shall make use of the dimensionless coordinate 
X = {Mloq/HY^'^ q, where M is the mass of the harmonic 
oscillator and q is the ordinary spatial coordinate with 
dimensions of length. In the Liouville space, we then 
have the pair of coordinates (x, x) to describe the system, 
which will transform as a hyperbolic rotation similar to 
A, cf. Eq. (35) below. Writing the A operators in terms 
of the position and momentum operators of the harmonic 
oscillator, we have 



x'±x' = e^^ (x ± x) , (36a) 
dx' ~^ dx' ^ \dx~^ dx] ^ (36b) 



where exp^ = \/b' jb. In Sec. IV A, we show that 
Eqs. (36) can also be generated by a symplectic trans- 
formation. 

In terms of the coordinate basis, the generator of the 
unitary transformation is 



G = i\ x—^-x— 
ox ox 



(37) 



A. Thermal symmetry in Fokker-Planck equation 

As a preparation to illustrate the existence of thermal 
symmetry in the classical Fokker-Planck equation [15], 
we first write the quantum MME in the coordinate basis 
before establishing its connection with the classical case. 

In the coordinate basis, the collision operator takes the 
form 



(31) 



K{x, x; b) = Kq{x, x) + Kd{x, x; b) 



(38) 
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where 
Kd{x, x; b) 





92 




dx^ 


\(- 


^ d\ 




dx J 




~{x 



, 2 ~2 
+ X — X 



(39) 







\dx 


dx) 



+ b- 



d_ d_ 

dx dx 



-{x- xf 



(40) 



arc written in a manifestly form invariant manner under 
the transformation in Eqs. (36). 

As an intermediate step, we introduce a new pair of 
coordinates, 



Q = \{x^x) 



(41) 



where the (r, Q) coordinate system is related to (a;, x) co- 
ordinate system by a counter clock- wise rotation of 45°. 
Eqs. (36) then imply that (r, Q) transform separately un- 
der the thermal symmetry as scale transformations 



Q' 



(42) 



where exp(6') = ^Jb' /b. 

To go into the Wigner representation, another pair of 
coordinates, {P,Q), where P is the Fourier conjugate 
variable of r, are used to describe the dynamics. Both 
coordinates transform with the same scaling factor under 
the thermal symmetry. 



Q' = e-^Q, 



P' = e-^P. 



(43) 



The collision operator in Wigner's representation is given 

by 



KiQ,P;b) = -.o{P^-Q^ 



d_ 
dQ 



Q 



d_ 
dP^ 



P 



92 



2 ydQ"^ 



d^ \ 
dP'^ J 



(44) 



As is evident from Eq. (44). the thermal parameter 6 
can be absorbed simultaneously into both coordinates, 
which is effectively a scaling transformation according to 
Eq. (43) . The existence of the second derivative terms in 
both the Q and P coordinates implies that a Gaussian 
wave packet would tend to a more uniform distribution 
along both the Q and P-axis as the temperature increases 
[23]. 

In Wigner's representation, the collision operator has 

the same form as the classical Fokker-Planck equation 
[15] expressed in the classical phase space variables, ex- 
cept with b replaced by its classical analog 



bci = ksT/huJo ■ 



(45) 



Therefore, we conclude that thermal symmetry exists in 
the classical system as well. We can then run the ar- 
gument in the reverse direction, starting from a given 
classical Fokker-Planck equation, we can transform the 
classical phase space variables corresponding to (P, Q) 
into the classical analog of the (a;, x) coordinate system. 
A thermal hyperbolic symmetry can then be defined on 
these coordinates as well. Indeed, the coordinates (x, x) 
in the classical system correspond to the Bargmann-Scgal 
representation [28], which is the classical analog of the 
coherent state representation in quantum theory [29]. 

Scale transformations of the form in Eq. (43) also apply 
in the coherent state representation of the MMEs [6], 
where the pair of coordinates z — P — iQ and its complex 
conjugate z* are used. A pair of "action-angle" variables 
(J, a) , where 



J=\{P^ + Q^), a = tan-i(Q/P), 



(46) 



can also be defined for the quantum system in parallel 
to the classical action-angle variables (for example, see 
[12]). Under the thermal symmetry, J then scales as 
exp(— 2^) J, while a is left unaffected. 



B. Thermal symmetry in Caldeira-Leggett 
equation 

In the previous discussion, we have focused on the ther- 
mal symmetry of the MME in Eq. (1). We will now ex- 
tend our consideration to other quantum MMEs. It is 
interesting to recognize that the Caldeira-Leggett (CL) 
equation [16] for the high temperature limit of a harmonic 
oscillator interacting with a thermal reservoir possesses 
this thermal symmetry as well. Indeed, the collision op- 
erator of the CL equation is given by 



Kci{x,x;bci) = Ko{x,x) - 71 (a; ~ ^) ( ^ ~ ^ 



libciix- 



(47) 



where bd in Eq. (45) is the high temperature limit of 
b. For Eq. (47), the second line (with coefficient bd) 
transforms under the thermal symmetry with an over- 
all multiplicative factor cxp(29) = b[,Jhcu whereas the 
other terms are invariant under the transformation. We 
therefore have the following transformation law for Kd, 



Kd{x,x;bd) Kd{x,x;b'd) ■ 



(48) 



and the CL equation is form invariant under the thermal 
symmetry. 



C. MMEs under more general conditions 

In order to determine the thermal symmetry for an 
MME, it is crucial that we are able to define for R{0) 
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one single angle that is independent of the frequencies 
of the field lo. This global property of the transforma- 
tion originates from the resonance effect expressed by the 
delta function J(a; — Wq) on the rhs of Eq. (6). Due to the 
delta function, h is restricted to the natural frequency of 
the reduced system wq- 

However, for MMEs derived under more general con- 
siderations (for finite temperature and without making 
the rotating- wave approximation) [17 19], the restriction 
of w to Wo does not occur in some of the coefficients. Let 
us consider the Markovian limit of the collision operator 
for the MMEs in Refs. [17-19]. They have the general 
form 



^ d d 

Kg{x, x; bk) = Ko{x, x) - 72(0; - x) [ — 



dx dx 

- 2^2b{x - xf + iT{x -x)[^ + ^ , (49) 



where 



Jo 



duj- 



(50) 



is a temperature dependent coefficient. The symbol P 
denotes the principal value of the integration and 
denotes the spectral density of the field. Unlike the coef- 
ficient 726, cf. Eq. (6), all frequency modes contribute to 
r due to the absence of a delta function in the integrand. 

The first 3 terms on the rhs of Eq. (49) have a structure 
similar to the CL equation and hence are form invariant 
under U{9). However, the F-term destroys the form in- 
variance of Kg for the following reason. To have the 
F-term remain form invariant, we require 



b'{oj) = e^%{oj) 
for all u), or equivalently, 

coth{haj/2kBT') coth{huJo/2kBT') 



coth{huj/2kBT) coth{hujo/2kBT) 



(51) 



(52) 



after using Eqs. (8), (7) and (25). This equation has no 
non-trivial solutions for T' > other than T' = T , for 
all u; > 0. 

We might think of introducing a transformation 
U[6{uj)\ with a frequency dependent angle, where 



6{u:) = \Hb'{uj)/b{u:)], 



(53) 



in order to restore the form invariance of Kg. However, 
since the field degrees of freedom have been traced out to 
obtain the MMEs, it is impossible to have all the U[6{uj)] 
(with different w) acting on the {x, x) coordinates simul- 
taneously. Therefore, this construction is not feasible. 

One might also consider modification of Eqs. (36) 
in order to incorporate a more general transformation 
that might account for the form invariance of the F- 
term. However, the requirements for the invariance of 
the canonical commutation relations in Eqs. (32) along 



with those for the form invariance of the first three terms 
on the rhs of Eq. (49) imder the modified transformation 
are too restrictive for this modification to be possible. 
For instance, if we modify the transformation coefficient 
for the derivative operator in Eq. (36b) from exp(±^) 
to exp(±(^), where 4> is an angle to be determined by 
the form invariance requirement for F, while we keep 
Eq. (36a) unchanged, we find that the second term on 
the rhs of Eq. (49) no longer remains form invariant. 
Furthermore, the modified operators no longer satisfy the 
canonical commutation relations in Eqs. (32). 

In conclusion, the thermal symmetry in Eq. (20) is 
applicable to MMEs of linear systems: (1) under the 
rotating- wave approximation with Eq. (1), (2) in the high 
temperature system with the CL equation (47) , as well as 
(3) in the classical Fokker-Planck equation. The symme- 
try does not apply to MMEs under more general consid- 
erations, for example without the rotating-wave approx- 
imation and for finite temperature. An example is the 
HPZ equation (49) [18]. For non-linear systems, gener- 
alization to the Bogoliubov transformations in Eq. (21) 
should be envisaged. 



IV. RELATIONS WITH OTHER 
TRANSFORMATIONS 

In this section we discuss the relation of the thermal 
symmetry to other forms of transformations utilized in 
closely related fields. 



A. Thermal symmetry as element of Sp(2n, K) 

Since the thermal symmetry preserves the commuta- 
tion relations of the A operators (13), it also preserves 
the canonical commutation relations for the position and 
momentum operators (32). We can write the commuta- 
tion relations in a more compact form by introducing a 
column vector 



X = 



(x\ 

X 

p 

\pj 



(54) 



The commutation relations in Eq. (32) can then be ex- 
pressed as 



[Xi , Xj 



*,J = 1,2,3,4^ 



(55) 



where is a non-singular, anti-symmetric matrix, 

in which I2 and O2 denote the 2x2 identity and null 
matrix respectively. 
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Eqs. (19) and their Hcrmitian conjugates can then be 
written as the action of an element of the symplectic 
group on X, 

X' = S-'^X, (57) 
where S G Sp(4, R) (with n = 2) is given by 



(58) 



We have previously introduced the rotation matrix R{9) 
in Eq. (21). The symplectic transformation S is gener- 
ated by 



J 



where 



I2 O2 



1 

1 



(59) 



(60) 



By construction, the symplectic group S preserves O 

SflS'^ = n , (61) 

as expected, where T denotes the transpose operation. 
^ The unitary operator U that corresponds to S acts on 
X by conjugation such that 



U{S)XU^{S) = x' = s-'^x 



(62) 



by Eq. (57). It can be expressed as U{S) = exp[i6'G( J)], 
where G( J) is the Hermitian generator of U that corre- 
sponds to J. It has the following expression [20], 



G{J) = ^Jf^njX = -{xp + l 



(63) 



which is identical to the generator of U{6) in Eq. (37) 

when expressed in the coordinate basis. The transforma- 
tion laws in Eqs. (36) then follow from Eq. (62). 

Notice that we started with a one dimensional system, 
in which the pair of canonical operators {x,p) is available 
for describing the dynamics. We are able to introduce 
the symplectic transformation Sp(4, ffi.) on the four com- 
ponents vector X only because the pair of tilde operators 
[x, p) are introduced when we go into the Liouville space. 



B. Relation of the thermal symmetry with SU(1, 1) 

There is a close relationship between U{9) and the 
SU(1, 1) group. Consider the Lie algebra su(l,l), which 
consists of 3 operators [21] 

M+ = At It ^ M_=AA, (64a) 
Mo = i(AtA-hItI+i). (64b) 



These operators satisfy the commutation relations 

[M_ , M+] = 2Mo , [Mo, M±] = M± . (65) 

An arbitrary SU(1, 1) transformation can be parameter- 
ized as 

with inverse U-^{9) = W{9) = U{-9). By choosing the 
angles as follows, 



7+ — —u , 



e.=9, 9o = Q, (67) 



Eq.(66) reduces to Eq. (17), or U{9) = exp{iG9), where 
the Hermitian generator 

G = i{M_-M+) = i{AA-A'<A'') (68) 

generates the Lorentz group S0(1,1) in the (a;, 5;)-plane 
[30]. 



C. Comparison with other form of the Bogoliubov 
transformation 



Let us compare the Bogoliubov (canonical) transfor- 
mation discussed in Ref. [22] with the Bogoliubov trans- 
formation in Eqs. (19). It was shown in Ref. [22] that 
MME with the following dissipative components 



where 



K^ = Kd + C3K3 + , 



Ksp = 2apa — aap — paa , 



(69) 



(70) 



can be recast in the standard form of the dissipative op- 
erator Kd in Eq. (1) by an SU(1,1) transformation on 

(a, at). 



(71) 



Here, u — cosh v and v — exp(z77) sinh i/ are coefficients 
related to ci , C2 and C3 . Even though this transformation 
reduces IC^ to the standard form K^, the uncoupled com- 
ponent of the dynamics Kq (which need not be Eq. (15)) 
may no longer retain its original form. 

In terms of the A operators, Eq. (71) becomes 



A = uA + vA^ , 



(72) 



which unlike the thermal hyperbolic rotation (20), does 
not mix the A with the tilde operators, A ot A^ . Hence 
this rotation can be reduced to a transformation in the 
Hilbert space. This shows that the SU(1,1) hyperbolic 
rotation symmetry of Eq. (71) is different from the ther- 
mal hyperbolic rotation symmetry considered in Section 
H. 
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V. OTHER ASPECTS OF THE THERMAL 
SYMMETRY 

A. Changes in the shape of Gaussian wave packet 
with temperature 

When the temperature of the reservoir changes, the 
profile of the density function p{x, x) changes in response 
[23]. If the temperature increases, the magnitude of the 
off-diagonal components of p decrease, a process called 
decoherence [9]. At the same time, the diagonal compo- 
nents of p stretch out to attain a more; uniform distribu- 
tion. The density function responds in the opposite way 
when the temperature decreases. With a Gaussian wave- 
packet as an example, we shall relate these behaviors of 
p through the thermal symmetry. 

Consider the cciuilibrium solution of the collision oper- 
ator C{b) in Eq. (28) in terms of the (r, Q) coordinates 
(41) [6, 11, 12, 23], 



1 



p{Q, r; h) = -= exp {-Q/2b - hr'' /2) 



(73) 



which is a Gaussian wave packet in both the r and Q 

coordinates. 

Wc now evaluate the averaged first and second moment 
of Q and r along the r = and Q = sections of p. 
Along the r = section of p, the average value of the 
first moment of Q vanishes, 



{Q)r 



J^^p{Q,0,b)dQ 



0, 



(74) 



where the notation (• • ■)r=o means that we are evaluating 
the average along the r = section of p. We also find 
that (r)Q=o = [35]. The second moments can also be 
evaluated readily to give 



((3^)r=C 



{r^)Q=, = l/b. (75) 



The dispersion of the wave packet (AQ, Ar) is defined as 

Ag = ^J{Q^)r=o - (0)^=0 . (76) 

with a corresponding expression for Ar defined for the 
(5 = section of p. Substituting Eqs. (74) and (75) into 
AQ and Ar, we obtain 



AQ = Vb, 



Ar = 1/Vb. 



(77) 



When the temperature increases, we have b' > b (or 
^ > 0), and the density function is transformed into 

p'iQ, r; 6) = p{Q', r'; b) = p{Q, r; b') . (78) 

The dispersion of the transformed Gaussian wave packet 
can be evaluated to give 



{AQY = Vb^ = e^AQ , 



'Ar, 



(79a) 
(79b) 



Eqs. (79) show that when the temperature increases 
{0 > 0), the Gaussian wave packet stretches away from 
the origin along the direction parallel to the Q-axis (r = 
section), whereas it contracts towards the origin along 
the direction parallel to the r-axis ((5 = section). The 
extent of the stretch and contraction is governed by the 
angle 0, or equivalently, the amount of change in the tem- 
perature. The density function behaves in the opposite 
manner when the temperature decreases. 



B. Disconnected regions in coordinate space 

When acting on the coordinate basis, the hyperbolic 
rotation R{0) preserves the bilinear form .x^ — x^. Each 
point on the {x, x)-plane belongs to a family of curves 
x"^ — = const generated by R{9) and these curves 
never cross. Consequently, the (a;, ^c) -plane can be di- 
vided into 3 disconnected regions, according to |a;| = \x\, 
\x\ > \x\ and \x\ < \x\. The fact that these regions 
are disconnected shows that the hyperbolic rotation does 
not mix quantum correlations (off-diagonal component 
(x|/9|x)) with probability (diagonal component {x\p\x)). 
For the one dimensional case considered in this paper, the 
|a;| > \x\ {\x\ < \x\) regions are further divided into 2 dis- 
connected pieces, depending on the value of sgn(a;) = ± 
(sgn(-r) = ±). For higher dimensions, besides an overall 
hyperbolic rotation dependent on the only available ther- 
mal parameter 6, there may also be ordinary rotations 
among the coordinates X = {xi,X2,---) orx= {xi,X2,---), 
without mixing between x and x. 

In analogy with special relativity, we can equally pa- 
rameterize the rotation matrix R in terms of a "velocity" - 
like parameter v. Indeed, let us define v by the relations 
[12, 14] 

cosh 61= l/Vl-v2, sinh6i = v/Vl-v2. (80) 
We find that v takes the familiar form 



V = tanh 6 ■ 



b + 



exp(-/zu;o/?) , 



(81) 



which is the relative probability of finding the frcqiiency 
mode Wo at thermal equilibrium. It asymptotically ap- 
proaches the "light" -speed v — > 1 (or — * oo), when 
T oo, or ft ^ in the classical limit. Hence, the 
classical limit as well as the high temperature limit is a 
singular limit similar to the speed of light limit in rela- 
tivistic theory. 



VI. CONCLUSION 

We have shown the existence of a thermal symmetry 
in the reduced dynamics of open quantum systems with 
non-unitary time evolution, for instance, in the MME 
with Eq. (1), in the Caldeira-Leggctt equation, as well 
as in the classical Fokker-Planck equation. However, for 
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systems considered under more general conditions, such 
as finite temperature and in the absence of the rotating- 
wave approximation, the thermal symmetry is no longer 
present. The Hu-Paz-Zhang equation is an example of 
the MMEs for such systems. 

The thermal symmetry is generated by a Bogoliubov 
transformation on the Liouville space of the reduced den- 
sity operator. This symmetry gives rise to degeneracies 
in the complex energy eigenvalues of the dissipative col- 
lision operator. The presence of this symmetry is due to 
the formulation of the dynamics on the level of the Liou- 
ville space; it does not exist on the level of the Hilbert 
space. 

As an important consequence of the symmetry, differ- 
ent thermal states of the system are connected, including 
absolute zero. Hence for systems observing the thermal 
symmetry, the effects of the thermal reservoir on the sys- 
tem can be inferred from the properties of the system 
at absolute zero, which does not contain thermal effects. 
From a different point of view, the effects of the sym- 
metry can be regarded as changing the value of Planck's 
constant for the case of fixed temperature. The symme- 
try then establishes a connection between the quantum 
and classical limit of the system. 

When represented in the coordinate basis, the thermal 
symmetry takes the form of a hyperbolic rotation on the 
dynamical variables of the reduced system. The angle of 
rotation depends on the amount of change in the tem- 
perature, and it governs the extent of stretch and con- 
traction of the density function along the diagonal and 
off-diagonal directions in the position coordinate basis, 



or the (x, a;)-plane. 

The thermal symmetry can be expressed in terms of a 
Bogoliubov transformation because the reduced system is 
a simple harmonic oscillator, i.e., a linear reduced system. 
Generalization of this symmetry to a nonlinear reduced 
system would require a generalization of the Bogoliubov 
transformation [31]. The thermal symmetry also enables 
the construction of a set of temperature dependent den- 
sity states of a coupled oscillator. We will present the 
discussion elsewhere. 
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